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Abstract 

We construct the equation of state (EOS) in a wide density range for neutron 
stars using the relativistic mean field theory. The properties of neutron star 
matter with both uniform and non-uniform distributions are studied consis- 
tently. The inclusion of hyperons considerably softens the EOS at high den- 
sities. The Thomas-Fermi approximation is used to describe the non-uniform 
matter, which is composed of a lattice of heavy nuclei. The phase transition 
from uniform matter to non-uniform matter occurs around 0.06 fm~ 3 , and 
the free neutrons drip out of nuclei at about 2.4 x 10~ 4 fm -3 . We apply the 
resulting EOS to investigate the neutron star properties such as maximum 
mass and composition of neutron stars. 
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I. INTRODUCTION 



The properties of neutron stars are mainly determined by the equation of state (EOS) 
of neutron star matter, which is charge neutral matter in /3-equilibrium. A comprehensive 
description of neutron stars should include not only the interior region but also the inner 
and outer crusts, therefore, the EOS for neutron stars is required to cover a wide density 
range. For the EOS at high densities, there are many efforts based on both non-relativistic 
and relativistic approaches, which discussed several possible mechanisms to soften the EOS 
at high densities, e.q., by hyperons, kaon condensates, or even quark phases [DHJ. When 
the density lower to 10 14 g/cm 3 , some heavy nuclei may be formed and matter becomes 
inhomogeneous. There are a few works based on non-relativistic models describing the EOS 
at low densities where the heavy nuclei exist |jJ-§]. Most studies of neutron stars are using 
the composite EOS, which is constructed by connecting the EOS at high densities to the 
one at low densities Even though the EOS at high densities are based on various 

relativistic many body theories, it has to be combined with some non-relativistic EOS at 
low densities. The differences in the models used in the different density ranges usually 
lead to some discontinuity and inconsistency in the composite EOS. Therefore, it is very 
interesting to construct the EOS in the whole density range within the same framework. 

In this paper, we provide a complete relativistic EOS for the studies of neutron stars, 
which is based on the relativistic mean field (RMF) theory. The RMF theory has been quite 
successfully and widely used for the description of nuclear matter and finite nuclei flOHl2 



We study the properties of dense matter with both uniform and non-uniform distributions in 
the RMF framework adopting the parameter set TM1, which is known to provide excellent 
properties of the ground states of heavy nuclei including unstable nuclei [fn|. The RMF 
theory with the TM1 parameter set was also shown to reproduce satisfactory agreement 
with experimental data in the studies of the nuclei with deformed configuration and the gi- 



ant resonances within the RPA formalism |Tj-|T6j. At high densities, hyperons may appear 
as new degrees of freedom through the weak interaction, the neutron star matter is then 
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composed of neutrons, protons, hyperons, electrons, and muons in /3-equilibrium. For the 
non-uniform matter at low densities, we perform the Thomas-Fermi calculation, in which 
the RMF results are taken as its input. The non- uniform matter is assumed to be composed 
of a lattice of spherical nuclei immersed in an electron gas with (or without) free neutrons 
dripping out of nuclei ]TT| , |r8|] . The optimal state at each density is determined by mini- 
mizing the energy density with respect to the independent parameters in the model. The 
phase transition from non-uniform matter to uniform matter takes place around 10 14 g/cm 3 . 
The same method (but without the inclusion of hyperons) has been used to work out the 
equation of state at finite temperature with various proton fractions for the use of supernova 



simulations [19,20 



This paper is arranged as follows. In Sec. II, we briefly describe the RMF theory 
and its parameters. In Sec. Ill, we explain the Thomas-Fermi approximation used for the 
description of non-uniform matter. The resulting EOS in the whole density range is shown 
and discussed in Sec. IV. We apply the relativistic EOS to study the constitution and 
structure of neutron stars in Sec. V. The conclusion is presented in Sec. VI. 

II. RELATIVISTIC MEAN FIELD THEORY 

We briefly explain the RMF theory used to describe the uniform matter. In the RMF 
theory, baryons interact via the exchange of mesons. The baryons considered in the present 
calculation include nucleons (n and p) and hyperons (A, E, S). The exchanged mesons 
consist of isoscalar scalar and vector mesons (a and u), isovector vector meson (p), and two 
strange mesons (a* and 0) which couple only to hyperons. The total Lagrangian density of 
neutron star matter, in the mean field approximation, can be written as 

£ = [ii^ - ( m B - 9aB<y - 9**bct*) - (g^B^ + g<pB<fr + 9 p bt 3 p) t° fpB 

B 
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where the sum on B is over all the charge states of the baryon octet (p, n, A, £ + , E°, E~, S°, 
and the sum on I is over the electrons and muons (e~ and fi~). The meson mean fields 
are denoted as a, u, p, cr*, and <fi. The inclusion of the non-linear a and u terms is essential 
to reproduce the feature of the relativistic Brueckner-Hartree-Fock theory and satisfactory 
properties of finite nuclei |T3 |. We adopt the TM1 parameter set for the meson- nucleon 



couplings and the self-coupling constants and some masses, which was determined in Ref. [JT3 
by reproducing the properties of finite nuclei in the wide mass range in the periodic table 
including neutron-rich nuclei. The RMF theory with the TM1 parameter set was also shown 
to reproduce satisfactory agreement with experimental data in the studies of the nuclei with 



deformed configuration and the giant resonances within the RPA formalism [|14]--|16|| . The 
hyperon masses are taken to be rri\ = 1116 MeV, mj = 1193 MeV, and mn = 1313 MeV, 
while the strange meson masses are m a * = 975 MeV and = 1020 MeV [Q. As for the 
hyperon couplings, we employ the following relations derived from the quark model 

\gaN = \gak = \qoy, = g^s, 

9 p n = \g P n = g P s, g P A = 0, (2) 

2(? ct *a = 2g CT *£ = g a *s = ^-g a N, g a *N = 0, 

2fiV,A = 2g^ s = = ^p-g^N, g^N = 0. 
In the RMF theory, the meson fields are treated as classical fields, and the field operators 
are replaced by their expectation values. The meson field equations in uniform matter are 
given by 

2 2 , 3 V- 2J b + 1 f kB rn* B 2 

m a a-g 2 a + g 3 a = ^g aB / k dk, (3) 



b 2n2 Jo y/k 2 + m% 

mlu + c 3 oo 3 = g^B (2J B + 1) 4/(6vr 2 ), (4) 

B 
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B 
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= (2Jb + 1) 4/(6vr 2 ), (7) 

where m* B = m B — g a B& — 9a*B^* is the effective mass of the baryon species B, and k B is 
its Fermi momentum. J B and I%b denote the spin and the isospin projection of baryon B. 

For neutron star matter with uniform distributions, the composition is determined by 
the requirements of charge neutrality and /^-equilibrium conditions. Considering the baryon 
octet and leptons included in the present calculation, the /^-equilibrium conditions, without 
trapped neutrinos, can be clearly expressed by 

Hp — Ht,+ — Hn~ He, 
Ha = Hz° — Hbp = Hn, 

HT,- — HS~ — Hn + He, 

Hfl = He, (8) 

where Hi is the chemical potential of species i. At zero temperature, the chemical potentials 
of baryon B and lepton I are given by 



Hb = yk 2 B + m* B 2 + g uB uj + g^ B <f> + g pB r 3 p, (9) 



H = y/tf + mf. (10) 
The charge neutrality condition has the following form: 

n p + n 2 + = n e + n n + n E - + n s -, (11) 

where ra, = kf/37T 2 is the number density of species i. Then, the total baryon density is 
tlb = n n + n p + rik + 7i£- + n s o + n s + + rig- + n s o. We solve the coupled equations (0)-(0), 
(H), and (JTTJ) self-consistently at a given baryon density The total energy density and 
pressure of the uniform matter are given by 

2J B + l f k B 
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III. THOMAS-FERMI APPROXIMATION 

In the low density range, where heavy nuclei exist, we perform the Thomas- Fermi calcu- 
lation based on the work done by Oyamatsu ||17||. In this approximation, the non- uniform 



matter can be modeled as a lattice of nuclei immersed in a vapor of neutrons and electrons. 
At lower density, there is no neutron dripping out of nuclei. We assume that each heavy 
spherical nucleus is located in the center of a charge-neutral cell consisting of a vapor of 
neutrons and electrons. The nuclei form a body-centered-cubic (BCC) lattice to minimize 
the Coulomb lattice energy. It is useful to introduce the Wigner-Seitz cell to simplify the 
energy of the unit cell. The Wigner-Seitz cell is a sphere whose volume is the same as the 
unit cell in the BCC lattice. 

We assume the nucleon distribution functions rij(r) (i = n for neutron, i = p for proton) 
in the Wigner-Seitz cell as 



rii (r) 



(nf - n° ut ) 



Ri 



+ n° u \ < r < Ri 

(14) 



n° u \ Ri<r < R, 



cell 



where r represents the distance from the center of the nucleus, and R ce ii is the radius of the 
Wigner-Seitz cell defined by the relation, V ce u = ^fRl e n = a 3 (a is the lattice constant). The 
parameters Ri and ti determine the boundary and the relative surface thickness of the heavy 
nucleus. R n and t n may be a little different from R p and t p due to the additional neutrons 
forming a neutron skin in the surface region. For neutron star matter at a given average 
density of baryons, % = J cell [n n (r) + n p (r)] d 3 r / V ce u, there are only seven independent 
parameters among the eight variables: a, n™,n°"*, R n ,t n ,n p n , R p ,t p . The optimal state is 
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determined by minimizing the average energy density, e = E ce u / V ce u, with respect to those 
independent parameters. 

The total energy per cell, E ce u, can be written as 

Eceii = E buU; + E s + E c + E e . (15) 

Here the bulk energy of baryons, Ebuik, is calculated by 

Ebuik = / Ermf ( n n (r) , n p (r) ) d 3 r, (16) 

J cell 

where Ermf is the energy density in the RMF theory as the functional of the neutron density 
n n and the proton density n p . As the input in the Thomas- Fermi calculation, Ermf at each 
radius r is calculated in the RMF theory for uniform matter with the corresponding densities 
n n and n p . The surface energy term E s due to the inhomogeneity of nucleon distribution is 
given by, 



E s = F Q \V(n n (r) + n p (r) ) | 2 d 3 r, (17) 

J cell 

where the parameter Fq = 70 MeV ■ fm 5 is determined by doing the Thomas- Fermi calcula- 
tions of finite nuclei as described in the appendix in Ref. ]I7 |. The Coulomb energy per cell 



Ec is calculated using the Wigner-Seitz approximation with an added correction term for 
the BCC lattice: 

E c =~ [ e [n p (r) - n e ] 0(r) d 3 r + AE C , (18) 

Z J cell 

where 0(r) stands for the electrostatic potential calculated in the Wigner-Seitz approxima- 
tion, AEc = CBcc{Ze) 2 /a is the correction term for the BCC lattice as given in Ref. |17 . 
n e is the electron number density of uniform electron gas, which can be determined by the 
charge neutrality condition as n e = Z /V ce u (Z denotes the proton number per cell). The 
last term in Eq. ([15]), E e , is the kinetic energy of electrons, which is given by 



where k e = (37r 2 n e ) 1//3 is the Fermi momentum of electrons. 

For each baryon density Ub, we minimize the average energy density e of non- uniform 
matter with respect to the independent parameters in the Thomas-Fermi approximation. At 
some higher densities, the heavy nuclei dissolve and the matter becomes homogeneous. We 
determine the density, at which the phase transition takes place, by comparing the energy 
density of non-uniform matter with the one of uniform matter. 



IV. PROPERTIES OF NEUTRON STAR MATTER 

In this section, we present the resulting EOS of neutron star matter in the density range 
from 10~ 7 to 1.2 fm -3 . At low densities where heavy nuclei exist, the non-uniform matter is 
described by the Thomas-Fermi approximation, in which the optimal state is determined by 
minimizing the average energy density with respect to its independent parameters. For the 
density below ~ 2.4 x 10~ 4 fm~ 3 , the nucleons form the optimal nuclei and those nuclei build 
a BCC lattice with uniform electron gas. It is found that the neutrons begin to drip out 
from nuclei at % ~ 2.4 x 10~ 4 fm~ 3 , then there is a neutron gas in addition to the electron 
gas. In Fig. [I] we show the neutron and the proton distributions along the straight line 
joining the centers of the nearest nuclei in the BCC lattice at the average baryon densities 
ub = 0.0001, 0.001, 0.01, 0.05 fm -3 . As the density increases, the optimal nuclei become 
closer and more neutron rich. At tib ~ 0.06 fm~ 3 , the nuclei dissolve and the optimal state 
is a uniform matter consisting of neutrons, protons, and electrons in /^-equilibrium. When 
the electron chemical potential exceeds the rest mass of the muon (at ub ~ 0.11 fm~ 3 ), it 
becomes energetically favorable to convert the electrons at the Fermi surface into muons, 
then the muons appear with the chemical equilibrium condition \x e = Hyperons appear 
at higher densities (ub ~ 0.27 fm~ 3 ). In Fig. |2| we show the fraction of species i, Y{ = rii/ns, 
as a function of the total baryon density tib- The composition of uniform neutron star matter 
is calculated by solving the coupled equations ©, and (|TT|). The threshold density 



for a hyperon species is determined not only by its charge and mass but also by the meson 



S 



mean fields, which are shown in Fig. |3| as functions of baryon density. In the present 
calculation, S~ is the first hyperon which appears at ~ 0.27 fm -3 , while A has almost 
the same threshold density (% si 0.29 fm -3 ). It is partly because the negative charge is 
much more favorable, even though S~ has somewhat larger mass compared with the mass 
of A. The other hyperons, S°, S + , H™, and S°, appear one by one at higher densities 
{jib ~ 0.57, 0.72, 0.84, 1.17 fm -3 ). The appearance of hyperons causes some decreases of 
the nucleon fractions. At high densities (jib ~ 0.7 fm -3 ), the A fraction is larger than the 
neutron fraction. We note that the hyperon threshold densities and fractions are sensitive 
to the hyperon couplings, and there are quite large uncertainties in the hyperon couplings. 
In this work, we adopt the hyperon couplings derived from the quark model. 

We display in Fig. |4] the pressure of neutron star matter as a function of energy density. 
The present EOS shown by the solid curve is compared with the EOS considering only the 
uniform matter phase (dotted curve), and it is found that the contribution from the non- 
uniform matter is quite large at low densities. The EOS without hyperons is also shown for 
comparison by dashed curve. The inclusion of hyperons considerably softens the EOS at 
high densities, because the conversion of nucleons to hyperons can relieve the Fermi pressure 
of the nucleons. In Fig. [5] we show the fraction of species i, Y>, in neutron star matter as a 
function of the average baryon density It is very interesting to see the phase transitions 
in the wide density range. At low densities, all nucleons exist inside nuclei, therefore the 
fraction of the nucleons in nuclei (dot-dashed curve) is equal to one. The decrease of the 
electron fraction (dotted curve), which is equal to the proton fraction due to the charge 
neutrality, implies that the optimal nucleus becomes more and more neutron rich as the 
density increases. Beyond the neutron drip density (jib ~ 2.4 x 10 -4 fm -3 ), there is a 
increasing fraction of the free neutrons outside nuclei (solid curve), and this causes a rapid 
decrease of the fraction of the nucleons in nuclei (dot-dashed curve). The phase transition 
from non-uniform matter to uniform matter occurs at ~ 0.06 fm -3 , where the heavy nuclei 
dissolve and the matter consists of neutrons, protons, and electrons in /3-equilibrium. We 
note that the neutron star matter is assumed to be at zero temperature, so there is no free 

9 



proton gas outside nuclei in the non-uniform matter phase. The muon fraction appears 
at tib ~ 0.11 fm~ 3 with the charge neutrality condition + Y e = Y p . At high densities 
(tib ~ 0.27 fm -3 ), the hyperon fractions appear, which have been shown more clearly in 
Fig. |. 

V. NEUTRON STAR STRUCTURE 

We calculate the neutron star properties by using the relativistic EOS. The neutron star 
masses as functions of central baryon density are displayed in Fig. ||. It is shown that 
the maximum mass of the neutron stars including hyperons is around 1.6M Q , while it is 
around 2.2M without hyperons. The neutron star mass is determined predominantly by 
the behavior of the EOS at high densities. The inclusion of hyperons considerably softens 
the EOS at high densities, therefore, results much smaller neutron star masses. The non- 
uniform matter, which exists in the crusts of neutron stars, has negligible contribution to the 
total neutron star mass, but it plays an important role in the description of the neutron star 
profile in the crustal region. In Fig. and |S|, we show the number density of the composition 
in the neutron stars with M = 1.6M Q and M = 1.2M respectively, as a function of radius. 
It is clear that the uniform matter containing the equilibrium mixture of nucleons, hyperons, 
and leptons exists in the internal region of the neutron star, while the non-uniform matter 
phase occurs only in the surface region. The neutron star with M = 1.6M Q has much thinner 
crusts as compared to the case of the neutron star with M = 1.2M . We show in Fig. ^|the 
mass-radius relations using the EOS with or without hyperons. It is found that the inclusion 
of hyperons only influences the neutron stars having large masses (M ~ 1.2M ). 

VI. CONCLUSION 

We have constructed the relativistic EOS of neutron star matter in the density range 
from 10" 7 to 1.2 far 3 . The non-uniform matter at low densities has been described by 
the Thomas-Fermi approximation, in which the nucleons form the optimal nuclei and those 
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nuclei build a BCC lattice. The uniform matter at high densities has been studied in 
the RMF theory. We adopted the RMF model with the TM1 parameter set, which was 
demonstrated to be successful in describing the properties of nuclear matter and finite nuclei 



including unstable nuclei ]nj, and its results were taken as the input in the Thomas- Fermi 
calculations. Hence we have worked out consistent calculations for uniform matter and non- 
uniform matter. The phase transition from non-uniform matter to uniform matter is found 
to take place at ~ 0.06 fm~ 3 . At high densities [ub ~ 0.27 fm~ 3 ), it is energetically 
favorable to convert some nucleons into hyperons via weak interactions. The inclusion of 
hyperons leads to a considerable softening of the EOS at high densities, since the conversion 
of nucleons to hyperons can relieve the Fermi pressure of the nucleons. We note that the 
contributions from hyperons are sensitive to the hyperon couplings, here we have adopted the 
hyperon couplings derived from the quark model. Presently, there exist large uncertainties in 
hyperon couplings. The hyperon couplings should be constrained by the experimental data 
of hypernuclei, but the experimental information is deficient for determining them. From 
the study of single A hypernuclei, the quark model values of A hyperon couplings usually 
predict overbinding of A single particle energies. It seems that the quark model values of A 
hyperon couplings lead to rather strong attraction. This might cause an earlier appearance 
of A hyperon. A detailed investigation of the dependence of the results on the hyperon 
couplings is deferred to future work. 

We have employed the present EOS to calculate the neutron star properties. With the 
appearance of hyperons, the maximum mass of neutron stars turned out to be 1.6M . It is 
found that the inclusion of hyperons results much smaller neutron star masses due to the 
softening of the EOS. The core of massive neutron stars is then composed of the equilibrium 
mixture of nucleons, hyperons, and leptons. The non-uniform matter exists only in the 
surface region, which forms quite thin crusts of neutron stars. The consideration of the 
non-uniform matter phase has negligible contribution to the neutron star mass, but it is 
essential to provide a realistic description of the neutron star structure. 

The present calculations have been performed within the framework of the relativistic 
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mean field approach, which is incapable to include pions explicitly. It will be possible and 
important to construct a complete EOS based on more microscopic theory such as the 
Dirac-Brueckner-Hartree-Fock approach. Especially, same approach should be employed 
in the treatment of both uniform matter and non-uniform matter. It is well known that 
the relativity plays an essential role in describing the nuclear saturation and the nuclear 
structure, it also brings some distinctive properties in the EOS comparing with the case in 
the non-relativistic framework. Therefore, it is very interesting and important to study the 
astrophysical phenomena such as neutron star properties using the relativistic EOS. 
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FIG. 1. The neutron distribution (solid curves) and the proton distribution (dashed curves) 
along the straight lines joining the centers of the nearest nuclei in the BCC lattice at the average 
baryon density n B = 0.0001, 0.001, 0.01, 0.05 fm~ 3 . 
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FIG. 2. The fraction of species i, Yi = rti/nB, as a function of the total baryon density ng. 
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FIG. 3. The meson mean fields as functions of baryon density. 
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FIG. 4. The pressure P versus energy density e for neutron star matter with the inclusion of 
hyperons (solid curve) and without hyperons (dashed curve). The EOS considering only uniform 
matter phase (dotted curve) is also shown for comparison. 
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FIG. 5. The fractions of the composition in neutron star matter as functions of baryon density. 
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FIG. 6. The neutron star masses as functions of central baryon density 



20 



FIG. 7. The number density of the composition in the neutron star with M = 1.6M as a 
function of radius r. 
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FIG. 8. Same as Fig. @ but for M = 1.2M ( 
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FIG. 9. The mass-radius relations for neutron stars. The solid curve shows the results with the 
inclusion of hyperons, while those without hyperons is plotted by the dashed curve for comparison. 
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